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Abstract 


A specific need for low da.ta rate coimnunication stiidy arises in ^ the siaeill probe 
communication link in deep space missions. The low available transmitter power aJid 
the large frequency uncertainty constrain the data rate to be low. An all- digital ^ 
communication receiver Is proposed' and its feasibility is estaolished. Although^ 
coherent systems should be used whenever practical, the noncoherent MFSK sy^^bem is 
more suitable for very low data rates. The effect of Rician fading on the performance 
of MFSK receiver is studied. Fading characteristics of the Venus channel are examine 
based on the exponential model and available experimental data on the Venus atmospliere , 
Because of the' requirement of high communication efficiency, three concatenated codes 
are evaluated and compared. The ra,pidly varying phase error at low data ra.te has 
great effects on the tracking loop behaviors which are examined by extensive computer 
study of the phase plane trajectories. Other toxucs discussed in this report include 
the spectrum of split-phase FSK and the coding/modulation selection for Pioneer -Venus 
communication systems . 
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I Introduction 


Chapter 

- surface of the near planets via landers as well 

The exploration of the ^■oriace oi 

flyby misBions toward the far planets will, for at least the next decade, he ^ ^ 

restricted dne to the weight limitations and the amonnt of electrical power a.arl- 

ahle. AS a result, the use of hiah power transmitters operating in conjunctron 

n Thus the power level of telemetry 
with a directional antenna may not he possi 1 . 

ev, will usually be below that necessary to perform coherent 
signals received on earth hill usual y 

aetection. The bit rate is low due to the power r.epulrements . 

a' possible solution to the problem, the use of M-ary noncoherent FSK as 
been considered for application in low-power space communications. Several forms 

„ , 

nf larpe freauency uncertainty caused largely hy dopp er 
detection in the presence of large ireauency 

^ ■ T.hiVitv The problems of time and frequency 

shifts and the oscillator instability. P 

, r' T '^,'rb'=‘rific multipath models, the 

synchronizations have been anaiyz^u L-i. -■ 

1 nrvofion svstem performance have been examined 

' effects of multipath on the communication system p 

'’■'’.a. m. • w.=a-i i— 

.SK i. « .... « ... »■. *“ 

ri<=-rform exactly the same as the optimum non- 
ideal case the FFT receiver should perform - . - y 

.1 -ir^rx vpite allows the use of 

TUPi c;v-< 5 tem operation at low data rape 
coherent MFSK receiver. The system op 

.ophisbicated Signal processsin,. The performance of bhe disibal system, however 
is limited by the effects of h^antization and finite word length. These problem, 
exam^ined in Chapter 3. Theoretical performance of the optimum noncoherent 

, n -M fading is considered in Cliapter 3. For 

MFSK receiver with and without mu 

.oael assumed, it is found that for a given average error 

pne 

probability, the required increase in the energy-to-noise density ratio to 
COB, bat the multipath fading is reasonably s,„all. particularly when the brt 
duration -IF filter bandwidth product is large. 
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One important application of the low data rate commimi cations is the 
Pioneer Venus coimaunication link, l^ie dense Venus atmosphere causes a severe 
fading on the communication link between the descending atmospheric entry probe ^ 
and the Eai’th. The most important atmospheric effect is due to the turbulence 
in the Venus atmosphere. Based on the latest Venus 7 and Venus 8 data', we may 
conclude now that fading effect is not as severe as early predictions based on 
the Venus 4 data. In Chapter 4 , we examne the following problems: the 

probability of different fa,ding levels relative to the free space as a fimctlon 
of altitude, the fading rate, the autocorrelation function and the spectrum of 
amplitude and phase variations, and the spectrum spreading due to turbiaence. 

An exponential turbulence model proposed by DeWolf [T] is used in the analysis. 

It is believed that the turbulent Venus atmosphere will not cause serious 
communication problem in the ^ forthcoming Pioneer Venus mission. 

TUe of +he M-.arv'- norceherent cba.rine]. is, however.^ very limited 

especially at low signal-to-noise ratio 18 ] . A.dditional coding implemented by 
concatenating the inner code and the outer code can provide low error rate 
without excessive decoding complexity. Performance of three di.fferent concatenated 
coding techniques is studied in Qiapter 5. The split-phase baseband format has : 
become increasingly important. In Chapter 6, the spectrum of the modulated signals 
with sp3.it-phase baseband is examined. The discussion of this chapter clarifies 
the ambiguities among all available results [9] . 

The performaice of coramandand telemetiy systems, useful in deep-space communi- 
cations, is frequently affected by the radio-frequency phase error which is intro-- 
duced at the point of reception by means of the carrier tracking loop. In low 
data rate communications, the phase error may vary rapidly over the duration of a 
signalling interval. If the input phase is varying according to a polynomial 
function of ti]ne, then the phase error can be reduced to zero in steady state by 
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using higher order loops. Exactly how the phase error is varying in time is not 
known. One way to model the phase variation is to assume that it he a sinusoidal 
function of time. Extensive computer analysis of the phase-plane trajectories ^ has 
shown that there axe threshold values for both the frequency and amplitude of such 
sinusoidal variation. This is the subject of Chapter T« 

In Chapter 8, we consider different candidates for coding and modulation and 
the interplex syst-em which provides a coherent system mthout relying on the .perfect 
reference signal. Finally in Chapter 9, the conclusions and recommendations for 
further study is presented. 
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Chapter II The Fast Fourier Transform Receiver 

^ ‘ Introduction 

The advent of large scale integration has suggested a new approach to signal 
processing problems. In the design of coinmunication receivers 5 the trend has been 
toward the increased use of digital circuitries. The speed of the present digital 
hardware, however, has prevented the data rate from being in the range of gegabits 
per second. For low data rate communication such as between 1 kilobits per second 
to 1 bit per second or lower, the digitial system is not only technically feasible, 
but also has potential cost reduction over the existing receiving systems . The 
low data rate communication system we consider here is designed primarily for deep 
space communications. It is /ilso suitable for other applications as long as the 
data rate is reasonably low. The digital, filters ertployed in the system can operate 
fast enough so that it may be time-shared among a number of channels. 

a'ouli cation^ , as the eommunicatiori between the earth and 

a probe in the ' atmosphere of Venus, the small transmitter power of the probe has 
constrained the data rate to be low. The large frequency uncertainty caused 
primarily by the oscillator instability and the' doppler effects requires a large 
IF filter bandwidth and thus the received signal -to -noise ratio is very small. 

To have an efficient communication, it is necessary to remove such frequency 
uncertainty. A small amount of frequency error which cannot be removed may be 
tracked by an automatic frequency control (AFG) loop. In this chapter we propose 
a fast Fourier transformer receiver where all components are implemented digitally 
except the frequency dovm conversion. To perform the mixing operation at 2.3 GHz, 
the typical operating fi'equency, by using digital circuitry is not feasible at the 
present time. Correlation operation in the kilo-hertz frequency range, however, 
can easily be performed digitally. 
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There are practical limitations on the digital systems. For example, the 
computer word length is finite and the A/D conversion error and the round off 
error are unavoidahle, Such limitations will "be examined in detail. 

II, System Co nfiguration and the FFT Operation 

Ihe propsed system is shown in Fig. 1. V/ith the exception of the frequency 
dovm conversion and the local oscillator, the receiving system can he iniplemented 
digitally. The input to the A/D converter can have a frequency of several hundred 
hertz to a few thousand hertz depending on signalling frequency and the uncompensated 
frequency ujicertainty due to oscillator instability and doppler frequency variations. 
The sampling rate should he at least t>ace of the largest frequency of the input 
signal. Experience has indicated that the desirable sampling rate be fo^ir or five 
times of the signal frequency. The digitized data may be recorded in magnetic tape 
to guard against the loss of data due to system failure such as loss of lock in time 
or frequency synchronization. Digital filtering of the digitized data reduces the 
noise in the received signal. This operation is optional if the signal-to~noise 
ratio is high but is essential if the signal-to~noise ratio is low. The discrete 
Fourier transform (DFT) of the data is performed by using the fast Fourier transform 
(FFT). The square of the magnitude of DFT is proportional to the power spectrum. 

The decision is based on the frequency of the largest spectral component. In the 
absence of noise, the receiver will always select the correct frequency. However, 
the magnitude of the measured spectral peak depends on the word timing error. In 
the presence of noise, time sync, loop provides the word tracking. The frequency 
sync, loop determines the frequency drift and provides an up-to-date estimate of 
the actual frequency. The local oscillator frequency can then be adjusted according 
to the frequency estimate. As the data rate is low, there is sufficient time for 
on-line, i.e. real-time, operation of the complete system to provide continuous 
frequency and time tracking and the signal -to-noise ratio improvement via digital 


filtering. 
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The fast Fourier spectral analysis can be very sensitive to the signal-to- 
noise ratio. Consider two FSK signals of frequencies 100 Hz and 200 Hz. Each bit 
of the received data contains one of the two signals plus an additive Gaussian 
noise. The bit duration is 1 sec. and the sampling rate is 1000 samples per second. 
The bit dui'ation ~ IF filter bandwidth product is thus equal to 100. For signal -to- 
noise ratios (S/W)^ = 0.1 and O.Ol, Figures 2 and 3'show the povrer spectra of a 
FSK signal plus noise. All spectral peaks are detected correctly in the signal- 
to-noise ratios considered. To avoid possible false spectral peaks which cause 
decision errors, it may be necessary to do some signal processing to increase the 
signal -to -noise ratio. By doubling the sampling fate we can use twice as many 
signal samples per spectrum. The acijS'Cent points in the power spectrum can he 
averaged to give [1], 


2K-1 




i=0 


2N-1 

I exp [- j i(n + l) ] 

i=l 


( 2 . 1 ) 


n = 0, 2, U,...,N 


Figure 4 shows the power spectra of the ''two— point” averages given hy Eq. (2.l) . 
Very little is changed in the spectral peal^;s but the noise is smoothed somewhat. 
Increasing the number of averaged points to, say, 4 and 10 may not have the 
desired improvement. This is illustrated in. Figure 5 for (S/N)^ - 0.1 and Figure^ 

6 for (S/k)^ = 0.01. - It appears that 4 is the maximum number of points that can 
be averaged to obtain any meaningful improvement from spectral averaging. Other 


considerations of the spectral averaging have been given by Winkelstein [10] , The 
use of digital bandpass filtering can also improve the power spectrum. 

The FFT receiver usually has more difficulty to detect the highest signal 
frequency component unless the sampling rate is much higher than such frequency 
([ 3 ]* The signal -to-noise ratios in Ref. 3 should be divided by a factor of 2). 

In spite of some practical problems with the FFT receiver as described above, the 
use of sophisticated signal processing a,t low data rate makes the FFT receiver a 
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very feasible system for deep space communications. Furtheimiore, the Fi’T hard- 
ware presently available hus a- speed several times faster than the use of iFT 
software. Such improvrnent in processing time from using digital devices Is 
particularly important for real-time applications. 

The probability of error of the fast Fourier transform receiver is the same 
as the optimum noncoherent MFSK. (multiple fi-equency shift lieyed) receiver. If 
sampling, quantization and finite word length effects are considered, the performance 
■will no longer be optimum. The sampling rate can always be chosen to be large 
enough to have negligible effects. Quantization and finite word length effects 
may seriously limit the use of the FIT operation. These problems will be discussed 
later , 

III. F requency Tracking in the I^iFSK Receivers 

A large ajiiount of frequency drifts that cause uncertainty may be removed by 
prior knowledge and prediction. However, other drifts will surely remain, and it 
is necessary to track these in any practical system. Several techniques of frequency 
tracking have been proposed, Goldstein's technique [1] obtains a frequency dis- 
criminator characteristic (S curve) by taking the difference, F^, of the two terms 
in Eq. (2.l) as an estimate of the current frequency error, 

V » T* "K* (2*2) 

n n+1 ‘ n 

whe re 

2H-1 - P 

'•n=lir .1 exp(- j ^ in)| 

1=0 

is the nth spectral line. Let T be the signal duration. The nonu.nal signal 
frequency here is (n + |-)/T, centered between two adjacent spectral lines. F^ is 
filtered and used to correct the local oscillator tuning as in any frequency -locked 
loon. The nrocedure is simple as it involves only FFT . One practical problem, 
however, is that necessarily equal to zero even if there is no frequency 
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drift and no noise ♦ Furthermore the method is suitable ior very small frequency 
drift which may not be the case in practice. 

The second teclmique due to Ferguson [2] is to use a weighted average of the 
k closest spectral components, where k is some small integer. If the received 
signal is norminally at frequency n/T*, then the estimate is 


F 

n 


C 


k/2 

I 


i=k~k/2 


a. 

1 


n+i 


(2.3) 


where C is a normalizing constant and the a^'s are a set of linearizing weighting 
coefficients. This technique of course is highly dependent on the choice of the 

a. ’s . 

The third technique due to Chadwick [11] is based on the spectral lines nearest 


to the observed frequency. Let be the spectral component of the observed 
frequency and and r ^ be the adjacent spectral lines with frequency l/T Hz 

larger and smaller, respectively, th'^n the observed frequency. The frequency 

estimate is 

r — r 
^ tl ~1 

n 2r T 


{2.k) 


which ten(3s to have less accuracy than Eq, (2, 3) and is useful only to very small 
frequency drifts . 

Ihe fouxth technique due to Simon [12] is an improvement of the third technique. 
Ttie frequency estimate proposed by Chadwick and used in a closed loop tracker 
suffers from the fact that it is biaaed. In fact Simon [12] was able to prove 
that a frequency estimator based on spectral estimates taken at integer multiples , 
of l/T cannot be unbiased. To be unbiased, the conditional mean of the frequency 
estimate must have the, linear variation. It turns out that all frequency estimators 
made up of spectral estimates taken at multiples of l/T have a conitional mean 
■ydiich has zero slope at the origin. Simon proposed an estimator which is derived 
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from spectral estimates taken at the adjacent 1/2T points » His estimate is 

_ ( 2 .! 

^n "fTr + Tj^) 

+ ''2 

Although this requires a slight Increase in spectral computation, the variation' of 
the conditional mean of this estimator is linear in the nel^ihorhood of the origin 
and thus presents a better estimator for constructing a closed loop frequency 
tracking algorithm. 


The fifth tedinique digitally implemients the automatic frequency control loop 
for frequency tracking (See Fig. l). Detail discussion of the digital AFC loop 

given in Ref. 2 where it is shown that the variance of the frequency estimate is 
smaller than that of Eq. time sync, loop can be implemented digitally 

in a similar manner as the AFC loop . 

IV. Effects of Quantization and Finite Word Length 

The A/D (analog to digital) conversion error and the round-off error of the 
digital AFC loop are descrihed in Ref. 3. It is shown that a commercially available 
1^1 bit register length is adequate to pro^/lde a signal-to-noise ratio of over UO dB 

in digital filtering. „ . . 

. The round-off noise has a more serious effect in the FFT operation. Let E^' 


be the quantization interval. For N-point FFT, the round-off noise variance with 


fixed- point arithmetic is given hy [13] 

which is proportional to N. The variance can be reduced by scaling, such as the 
multiplication factor of 1/2, at each stage of the FFT operation. Scaling, of 
course, requires a slight increase in the di.gifal hardware. With scaling, the 
variance becomes 

? 

E" 

. ^ 12 ^ 


( 2 . 6 ) 


(.2.7) 
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2 

but the noise -to-signal ratio is 5WE which is still proportional to N. In low 
data rate communication, a typical value for N is 1024. To achieve a signal— to- 
uoise ratio of 40 dB, the required register length is 13 hits. If the siganl 
dynamic range is also considered, much longer register length is needed. Thus 
the number of points that can effectively be performed by a given FFT hardware or 
software is limited. If the floating-point arithmetic is used, the mean square 


( 2 . 8 ) 


error (round-off noise variance) is upper bounded by [l4] 

where K = 2^^ is a power of 2, For K = 1024, the variance given by Eq, (2.8) is ■ 
greater than that of Eq. (2.7) but smaller than that of Eq. (2,6), 

V . Concluding Eemarks 

V/ith the present technology in digital devices, we have shown that the proposed 
all-digital low data rate conmrunication system which employs a fast Fourier trans- 


fo 


iiii recej- ver j_S 


completely feasxble from eoth perioriiiancc and econoiiij 


The effects of the quantization and finite word length on the digital system per- 
formance are normally very significant. It is shown in this chapter, however, .that 


by properly designing the digital, filters and iinplementing the FFT operation, such 


effects can be minimized. A commercially available l4 bit register length can 
provide a signal -to -noise ratio of over '40 dB in both digital filtering and FFT. 
With continued improvement in the digital hardwares, there is every reason to 
believe that the proposed all-digital system provides cost-reduction over existing 


receiver systems. 
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Chapter III Performance of Wideband Noncoherent MFSK System with Multipath Fading 
I , Introduction 

Reports from the Mariner V S-band measurements [15], [l6] and the Soviet 
Space Probe Venera U, 5 and 6 have both indicated the severe fading of radio 
signals at the turbulent Venus atmosphere. Refraction of the radio beam by the 
charged particles of the upper atmosphere and the gases that constitute the lower 
atmosphere produce changes in frequency, phase and amplitude of the signal received 
at the deep— space tracking stations on the Earth. In addition, the amplitude of 
this signal was affected both by defocusing and absorption in the lower atmosphere. 

It is necessary to consider the performance of the noncoherent coded system under 
various fading conditions, Glenn [IT], Schuman [l8], and Chadwick [5] have 
examined the effects of multipath fading on the low data rate communications . 

Their results, however, are inconsistent and limited to the wideband binary FSK. 

The Gaussian assumption made by Glenn for very large IF filter bandT^/idth (b^Hz) 
and bit-duration (T^) product may not be valid for intermediate values of b T , 
say 1 < b T^ < 10. In this chapter the exact error probability of the wideband 
noncoherent MFSK receiver is derived from using the well known "Rician" channel 
model [19]. The received signal consists of the spectilar and the random scatter 
(diffuse) components. The results correspond to the situation that the reflected 
siganl has a much larger bandwidth then the direct signal and that the time delay 
between the direct and the reflected signals is much less than the bit duration. 

This represents the most in^ortant kind of multipath fading. Other multipath 
fading conditions that may limit the performance more but have a smaller probability 
of occurence are also considered. They are: (l) the reflected signal amplitude is 
assumed to be const ant over a bit period but varies randomly according to the 
Rayleigh distribution from bit to bit, and (2), the delay of the reflected signal 
exceeds one bit period. The reflected signal may then be considered as part of 


the additive noise. 
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II . Error Probability Computation 

Consider a noncoherent MFSK receiver which has 2^ parallel channels where n 
is the number oT information bits per word. In the absence of multipath fading, 
the probability density of the output of the kth channel which is mateched to the 
kth signal input is [20] 


p(Xj^) = 


2N b 
o o b 


exp - C 




+ 2nE 


/2nE 




2N b 
o o b 


^ ^ N b 

o o b 




> 0 


C3,l) 


= 0 otherwise 


where E is the signaH. energy per bit. For and with multipath fading. 


Lindsey [21] has shown that 

\ + p 

1 + B 




XV 


> \ 


> 0 


(3.2) 


(1 + 3) 

= 0 otherwise 

and the probability density of the output from other channels which are not 

•s 

matched to the kth signal is 

>0 j j k 

= 0 otherwise 

The parameters p and 3 in Eq . (3.2) are, respectively, n times the energy -to-noise 
ratio per bit of transmitted information produced by the specular component and 
n times the energy -to -noise ratio per bit of transmitted information produced by 


p(y^) = p(y) = (-y); y 


the scatter component, i.e. 

, a^E ^ ^ r ^ 

p = n( ) ; 3 = n( -r; — ) 


N 


(3.h) 

2 


o o 

where a is a factor proportional to the strength of the specular con^onent and cr 

is the variance of the scatter component. By comparing Eqs . (3.1) and (3.2), with 

b = 1, we note thht in going from "no fading’' to "with fading", the following 
o b 

change is made: 

E E ^ 1 _ 

N 1 + 3’^ \ l+ 3 
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By using the transformation x ~ — and multiplying ^ 

o 

Eq. ( 3 . 2 ) hecomes r -i 




= 0 otherwise 


; X > 0 


(3.5) 


At the output of the kth channel, the sample average of M - samples is denoted 

as u. The characteristic function of u ia ^ , 


(3.6) 


X ,M “ N 1 + e ■ 1 - i(l + K)t'^ 

{1 - rCl + pit] o J 

Similarly the sample average v of M samples in other channels has the characteristic 


function 


4>(t) = ■ '-'3f 

(1 - it)^ 


(3.7) 


The probability densities p(u) and p(v) can be determined from Eqs . (3.6) and (3.7) 
respectively. The exact error probability is 



2^-1 

1 p ( V . ) dv . 

J J . J 

Pj,(n) = 1 - I p(u)du 

jT 

j=l 

0 

0 


^ r (2 ^ ^ ul r 

■S^-i /n 1+p ^ 

f nE : a N ° ' ' (M-l) ' ' 1 

= 1 - f No 1 +g^ JnEoi jM-l ' ^ 


^ J 


ru(i+p) 


M-l -V^ 

V e dv 


Equation (3.8) is consistent with Lindsey's result [21] for M = 1 and Chen's result 
[20] for no fading case (p " 0, a = 1), 


I (2 ^) 

_ / ^ . / -(u+;: 7 “ ) M-1 ^ o 

P„(n) = 1 - e No u 

^ , nEu xM-1 

' o 


\ 

1 M-1 -V, 

wTT i ® 

0 J 


2"-l 


1 -V M“1 , T 

CM‘-ll! j “ ' “■ * 

o 


M-1 M-k-1 

-u Y _u 





A listing of computer programs to calculate Eq. (3.8) is given as Appendix of this 
chapter. Equation (3.8) is also similar to an error probability expression derived 
by Lindsey ([22], Eqs . (37) and (33)). Accurate computation of the integral in 
the Lindsey's expression has recently been made by Adams [23], which is adapted 
to the computation* of Eq, (3.8). For M = 5, 10, and 100 with N = and 64, the 
error probabilities are plotted in Figs. 7, 8, 9 respectively, versus the total 
average received signal energy per bit to noise density ratio, 

(1 + Y^)&/n = (a^ + 2o^)E/K^ (3,10) 

2 5 2 2 

for Y = 0,1,10, and 10 (^o). For y = “=, (l + y ) B/n = E/N^. It is interesting 

to note that for communication channels which are largely scatter in nature, i.e. 

Y <2, the error performance improves as M increases for large • Also the 

perfomance degradation due to multipath fading decreases as M becomes large. This 

is probably due to the fact that the MFSK receiver is a noncoherent energy detector 

s 

and the reflected signal tends to increase the total signal energy. 

III. Other Multipath Fading Conditions 

The above results are based on two, assumptions: (l) the reflected signal is 

purely random and the direct signal contributes to the nonrandom (specular) component 
of the received data, and (2) the reflected signal has a negligible time delay from 
the direct signal. The second assumption is clearly justified in view of the vast 
communication distance involved. The first assumption holds when the altitude of 
the space probe above the Venus surface is large enough, that the Venus atmosphere 
has little effect on the direct signal. 


*The author would like to thank Dr. Lindsey for calling his attention to the 
computer program prepared by W.B. Adams, and Mr. Adams for using his computer 
program [23] . 
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If the second assumption is not satisfied, the reflected signal may not he 
received at the same bit duration as the direct signal. Consider the worst case 
that the reflected signal provides no infoimiation and can be included in the 
additive noise. Then the .average received energyr^to-noise density ratio per 
information bit is 

g^E ^ ■, (3.11)- 

2a^E + N 1 + 

o P 

That is the total received energy -to-noise density ratio per bit is reduced by a 


factor 

(1 + Y^)(l g/n) > 1; 0 (3.12) 

The error probability given by Eq. (3.8) is still valid if the required increase 


given by Eq, (3.12) is taken into, account in plotting the error probability. For 
= 0, then the receiver does not receive any useful information and the 
probability of error is 1 - l/N. As 3 increases, Eq. (3.12) indicates that the 


required m or ease 


energy — oO— nois( 




1 ■i'OpvjYY. -Pinn or 


of 3 . 

Nexb we consider the case that the first assumption is not valid. The re— 

% 

fleeted signal is now assumed to be constant over a bit period but varies randomly 

according to the Rayleigh distribution from bit to bit. Furthermore, the time 

delay of the reflected signal is assumed to be negligible. The error probability 

can be derived as follows. First obtain the error probability without fading by 

setting 3 = 0 and a = 1 in Eq. (3.8). Consider now S = E/T^ as a random variable 

with Rician distribution, 

p(S) - exp i 

2c 2c 

"s. 

where S is the power of the direct signal component. The error probability without 
b 

fading is then multiplied by Eq. (3.13) and integrated with respect to S from 0 to 


(S + S ) — 

C 
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In conclusion, it is remarked that multipath fading i^ an important 
consideration in the forthcoming deep space missions . An accurate eva,l-uataon 
of the multipath fading effect requires a realistic channel model which has to 
be derived from the experimental results. As severe fading occurs only for a 
short period of the entire mission, the percentile error probability must also 
be conputed. • 


-IT- 


Chapter IV Atmospheric Effects on Pioneer Venus Communication Links 
I, Introduction 

Based on the communication signals received from Venera U, 5, ^ n^d 7? there 
is definitely observed fading which may he caused hy turbulence in the Venus 
atmosphere. In order to determine the possible effects that turbulence could have 
on an S-band communication system, theoretical studies have been made [T, 22, 23, 

2t, 25, 26 ] to compare with the limited e:xperimental data available. Assumption 
was made in these studies that the observed effects were indeed due to turbulence. 

It was recognized, however, that other factors such as absorption, refraction and 
defoc using, and the une^jplained motion of the probe, may also contribute to the 
observed fading. But since turbulence appears to play the dominant role, it is 
necessary to see what type of S-band effects mi^t produce, and the implications 
of these effects on the S-band communication systems. 

Venera-T data indicates that there is no fading at 1 GHz with radio waves 
propagating vertically through the entire thickness of the Venus atmosphere. 

Fading is still severe as the space probe traversed the planet's atmosphere. 

These conclusions are not consistent with the earlier results from Venera-^ data. 
Venera-T data should be more reliable, however, as the probes did land on the 
Venus surface. Additional reports on turbiiience effects and Venera-T data are 
available in Refs. 2T-31. Reference 2T is particularly informative on Venera-T 
data. A summary of U.S. efforts on Venus study is given by Ref. 32. 

Along with the limited experimental data, theoretical studies are useful to 
draw meaningful conclusions of the turbulence effects on radio wave propagation 
through the Venus atmosphere. The book by Tatarski [33] provides a useful back- 
ground for the study. The theoretical problem areas to be studied are as follows: 
1. The probability of different fading levels relative to free space as a function 
of altitude. 
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The fading rate, i,e. the number of fades per second. 

3. The autocorrelation function of the amplitude fluctuations about the mean 
value , 

The autocorrelation function of the phase fluctmtions about the mean value. 

5. The phase and frequency spectra. 

6. The spectrum of the frequency fluctuation about the carrier frequency, i.e. 
the spectrum spreading due to turbulence. 


II . Fading Levels 

Fig. 10 shovs the geometry assumed for the propagation problem. The procedure 
for calculating fading levels is as follows. 


The fading of Venera— 4 as a function of altitude can be approximated by [26], 
a^(z, 0.94 GHz, 6 = 0'^) = 1.175 exp(- O.O 785 2) . (4.l) 

over the 20-40 Km altitude range, where 



(4.2) 


A is the signal an 5 >litude, z is the hei^t above the Venus surface, and <*> denotes 


the time average. These' .measurements were taken at a frequency f = 0,94 GHz and at 
a sub-earth point, i.e. e = 0°. Note that the scale height as defined by Eq. (4.1) 
is h = 1 / 0.0785 = 12,7 Km. When the frequency and 0 are changed, scaling laws 


must be \ised to calculate a^. Under reasonable assumptions, the following I26J 

2 

are two possible cases for scaling o . 


Case (a): 



(4,3) 


Case (b): L 

o 

2 

• 



(cos0) 


(It.l*} 
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(!*.7) 


where is the outer scale of turbulence, k = 2 WA . and A' is the wavelength. 

Using the scaling laws (4.3) and (4.4) and Eq. (4.1), we find. 

Case (a) o/z, 2.3 GHz, 6 = 60°) = 4.07 exp(- 0.0785 s) 

Case (b) o/z, 2.3 GHz, 6 = 60°) = 3.74 exp(- 0.0785 z) (**• 

, Given a^, the fading level distributions can be calculated assuming that the 
probability distribution of the amplitude scintillations is Rice-Nakagami (a 
constant vector plus a Rayleigh distributed vector). Based on DeWolfs for^^la 

mb, • T TTT-i-v. field E(s) becomes a constant-- 

[22], we let . The signal with electric iieiu \ 

plus-Rayleigh distributed vector, i.e. 

E(s) = E^(s) [exp(-a^) + <5B] ^ 

where 6 B is Rayleigh distributed, with zero mean and <| 6 b|^ = 1 - 2o^)- 

por d or a calculated in (4.5) or (4.6), we can compute the Horton parameter K 
A e 

defined by DeWolf [ 22 ] as 

2 

K = 10 logj^Q[exp(2o^) - 1] 
which is the ratio of powers in the two terms of Eq. (4.7) 

Using the formulas in Norton's paper [34], the prohabilxty of different fadi g 

levels may be found. A computer program which was originally prepared by 
Dr. Strohbehn was adapted here for the computation of the fading levels for various 
probabilities as a Action of altitude. A listing of the computer programs for 
fading level computation is given as the Appendix to this chapter. Because of the 
uncertainty on L^, both Eqs. (4.5) and (4.6) were used for scaling. For carrier 
frequency 2.3 GHz and 0 = 60°, the two scaling laws provide essentially the same 
results. For probabilities less than 0.05, the fading level vs the altitude plots 

Shown in Fig. U have dips at altitudes between 20 and 25 km for both scaling laws. 

. _ . . i . .« 4 - 4 -'u„ 4 , -Pfvr /T 0 n CTabin g "ths 

The problem [3^ is essentially due uo une .aiio 0440 .. 44 


(U,6) 
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fading levels are polynomial approximations in the parameter ^2 

computer printout). Since only a few terms are kept in the approximations, the 
results are very poor for k values close to unity. It is noted that the dips 
all occur between 20 and 25 kra where the k value is going through unity. A 
simple curve fitting through this trobulesome region still provides accurate 
results. Fig, 12 is the fading level plot for various probabilities. .The tran- 
sition from p = 0.15 to p == 0.20 is not a smooth one; so two separate plots are 
provided in Fig. 12, Other fading level plots are given by Chen [S6]and 
Strohbehn [26 ] . 

Ill, Fading Rates 

Given that there may be strong fading, the next parameter of interest is the 
fading rates. In order to get a feel for the magnitude of this parameter, it was 
approached from several different directions . The first two estimates were based 
on the observed fading rate of 1 fade/sec. as measured by Venera-H and reported by 
DeWolf [22J. In the first estimate, it is assumed that the fading rate will scaJ.e 


as 


f ) - 


C 


f (f ) 

C. 1 


(4.9) 


where f is the fading rate and f and f^ are the carrier frequencies. This 
A ^2 

formula assumes the fading is dominated by local winds. It gives a value of about 
1,5 fades /sec, for a carrier frequency of 2.3 GHz. The second estimate assumes 
that the motion of interest is the velocity of the probe, and the scaling law 


becomes 


Again at 2.3 GHz, “ 1*5 



v^/vi where v^/v^ 


(4.10) 

is the ratio 'Of the velocities of the 


new probe to the Venera probe. 



- 21 - 


The third, estimate was made completely independently *of the Venera-4 data. 

It is known from turbulence theory that the critical size of an inhoinogeneity to 
produce fading is roughly , 

^ _ /ttAiT (4.11) 

where L is the effective path length in the turbulent medium. The fading rate then 
is 


V 

f = .£. 

A 

where v is the velocity of the probe perpendicular to the line-of-si^t of the 
P 

propagation path. For the case of interest here, 

0 . 866 V 


(4.12) 


■ = 


(It .13) 




ttXL 


For X = 15 cm, we get f^ = 3.6 v/vfT. If the probe is moving on the order of 
30 m/sec in the lower part of the atmosphere, then we get ~ 0.54 to 1.5 fades/sec. 

as the effective length, L, varies between 40 Km and 5 Km. 

Ihe above calculations give us a rough idea for the e:xpected fading rates. 

A more sophisticated model can be used to predict the fadi.ng rate as a f vine t ion 
of altitude by knowing the velocity of the .probe and' using sxi exponential model 


for the turbulent atmosphere. 

For Venera-T, the vertical velocity before impact is IT m/sec. The fading 
rate from Eq.. (4.12) is O .985 fade/ sec. It may be concluded from the above 
discussion that the maximum fading rate is 1.5 fades /sec. 


XV. Concluding Remarks 

The autocorrelation functions of the amplitude and phase fluctviations , and 
the phase and frequency spectra and spectrum spreading are all described in detail 
in Refs, 26 and 36 . It was shown that neither the pahse fluctuations nor the 
frequency spreading would be a major problem in system considerations. The very 
limited preliminary report available on the ”Venus- 8 ’* data indicates that the 
data from **Venus— 8 ** are similar to that from Venus— T [35] * 
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It is reasonable to conclude that based on the Venera data and the 
exponential turbulence model proposed by DeVolf, there will always be strong 
fading at low altitude (below 20 Ion), particularly near the Venus surface. 
Fading is least significant if 6 = 0°. Further study on the subject matter of 
this chapter is much needed, however, to assure the success of the forthcoming 


Pioneer-Venus mission. 



-23- 


Chapter V Concatenated Coding for Low Data Rate .Comniimi cations 

I . Introduction 

In deep space communications with distant planets, the data rate as well as 
the operating signal-to-noise ratio may he very low. To maintain the error rate 
also at a very low level, it is necessary to use a sophisticated coding system 
(longer code) without excessive decoding complexity. The concatenated coding has 
been shown to meet such requirements in that the error rate decreases exponentially 
with the overall length of the code while the decoder complexity increases only 
algebraically. Three methods of concatenating an inner code with an outer code 
are considered. Performance comparison of the three concatenated codes is made. 

It is shown that the concatenated code with inner code a convolutional code and 
outer code a Reed-Solomon code performs the best among the three. Refs. 3T-^6 
contain most of the available informations on the Concatenated code. Performance 
of the concatenated codes in the presence of multipath fading has not been con- 
sidered. This chapter is based on Ref. 

II. Concatenated Codes 

The three concatenated' codes considered are: 

Code I. inner code a bi-orthogonal code, outer code a generalized Hamming code 
or Reed-Solomon (R-S) code. 

Code II. inner code a convolutional code, outer code a block orthogonal code (?4FSK) 
Code III. inner code a convolutional- code, outer code a R-S code. 

Other methods of concatenation are possible. For example, two k = 6 by v = 2 
convolutional codes can be concatenated to five a 9 by 4 code . Here k is the 
constraint length and l/v is the rate of such a code. Each convolutional code can 
be decoded by using Viterbi’s decoding scheme. Erickson [i+l] has shown, however, 
that the concatenation of Viterbi decoders does not appear to be useful in the 
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present context of tlie planetary program. He conjectured that the most 
appropriate outer code, in any concatenation scheme involving a Viterbi algorithm 
inner decoder, vas a high-rate algebraic block code. 

Consider first a generalized Hamming code as the outer code [38]. The code 
has n elements including k information (data) elements and m = n - k check elements 
Each data element is a six-bit bi -orthogonal code word. The receiver performs 
both error detection and correction. The generalized Hamming code which has a 
Hamming distance of three is a specific case of R-S code. For the R-S code, the 
minimum distance d, between two code words is related to the number of check 
elements, m, by d = m + 1 . The maximum number of correctable elements for each 
code word, t, is equal to m/2. The probability of the bit error after both error 


de'tection and correction is 

- ^ (5.1 

c - Ag + A,, 

where = probability of bit error at the detector outpun , - prooabilxu^ Ox 

bit error at the corrector output, and and are the normalized data quantities 
at the detector and the corrector outputs respectively [38] . 

Without restricting to the Hamming distance of 3, Forney [3T] Simpson [40] 
have considered the bi-orthogonal inner code and the R-S outer code. A typical 
concatenated coding system is shown in Fig. 13 where the R-S code can correct up 


to 2 errors. The main difference among the three reports [371, [38], [40] , is in 
the decoding method. Forney considers both the maximum likelihood decoding and 
the generalized minimum distance decoding of the R-S code. Both Miller and Simpson 
use the algebraic decoding for the R-S code although their error probability 
expressions are inconsistent . The digit error given by Simpson is 

? /n\ ^ '^-i 

^ ii/ 
i=n-t^ / 


n-t-1 

Pj.Ce) I 

^ i=t+l 


f+. + -i^ AnN^i _ ^n-i 

X — f I X I_1 II 


n 




F U - F ; 
e e 




p 
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where P is the probability of error of the bi-^rthogonal code word. Available 

C . ' • • 

results are sho^ in Fig. l4 with the error probability plotted as a function of 
signal energy per bit to noise density ratio of the inner code. Curve 1 has a 
2^ - symbol code and a. bi -orthogonal code of rate l/l6 [3T]. Curve 2 has a 
(l8, 12) R-S outer code (12 information elements out of a total of l8 elements) 
and a 6 bit bi-orthogonal inner code [38]. Curve 3 has a ( 63 , 49) R-S outer code 

and a 6 bit bi-orthogonal inner code [40]. Although the code lengths are not the 

same, it is clear that Curve 1 is better than Curve 2 which is better than Curve 3. 

■ The block diagram of the concatenated system for Code II is shown in Fig. 15. 
Although it is possible to determine the upper bound of the error probability, the 
bound may be too loose to be useful. Some computer simulation results were reported 
by Richardson et.al. [44]. Let M be the bit-duration - IF filter bandwidth product. 
Figure 16 is a plot of the error probabilities for v = 3, k = 6 inner code and 
V = 5 , k = 8 inner code with M = 2 and 10 along with the performance of the coherent 
systems without concatenation. The degradation in performance M - 2 to M = 10 
is approximately 1.1 to 1.4 dB in the signal energy per bit to noise density 

ratio. It is noted that the performance improvement over the wideband noncoherent 
MFSK system [20] is very significant for both M = 2 and 10. 

Code III was considered by Odenwalder [39]. A block diagram of the system is 
shown in Fig. IT. For any inner code of contraint length k bits, there exists a 
R-S outer code of block length 2^-1 and with 2^ - d information bits- If is 
the error probability of the convolutional inner code, then the word error probab- 
ility with the R-S outer code is 


P 

w 


(e) 



P^ (1 - P 
e e 


(5.3) 


where d ^ 2 t + 1 and 


is the closest integer to d/2. Althou^ analytical expres- 

— 1“ j 

Sion of error bound of Code III is available for very limited cases [46], the com- 
pute.r simulation result Is more useful. The best computer simulation result of 
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g 

(l/3, 8) convolutional Inner code and 2 - symbol R-S outer code is shovm in 

Fig. l8. Here y ~ 3, = 8. Also plotted in the same figure are the error probab- 

ilities of a 5-bit bi -orthogonal code, an (l/2, 8) convolutional code vith the 
Viterbi’s maximum likelihood decoding (from Odenwalder [39], Chapter 5), and a 
concatenated code (31, 25, 5) with (31, 25) R-S outer code and 5 bit bi -orthogonal 
inner code. The improvement from bi -orthogonal only to bi -orthogonal /R-S code is 
added to the (l/2, 8) convolutional code to give an estimate of the error prohabillty 
of a concatenated code with (l/2, 8) convolutional inner code and (31, 25) R— S code. 
The error probability of Code III is the best among all concatenated codes . The 
required interleave-buffer, however, essentially increases the word length, or 
reduces the effective signal -to -noise ratio. 

Performance curves indicate that the concatenated codes greatly .improve the 
peformance over the uncatenated codes. The decoding complexity increases. from 

s 

Code I to Code II to Code III while _the performance improves in the reverse order. 

The decoding complexity depends manily on the inner code used. For the low 
available signEil-to-noise ratio, Code III is definitely the best choice to keep 
the error rate also at a low level . 
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Chapter VI On the Spectrum of Split-Phase FSK 


I. Introduction 

It is now well known that the use of a split-phase baseband signal provide^ 
better DC response as compared with the use of a non return-to-zero (NRZl baseband 
signal. The split-phase signal, however, req^uires a larger transnaasion bandwidth. 
The discrete frequency components that always exist in the spectrum of split-phase 
FSK facilitate acquisition and tracking of bit clock. Spectral analysis: of the 
split -phase FSK has been reported by Hartmann [^8], Shehadeh I^9J and Chen I50J ,. 
Such analysis is needed to determine the bandwidth of the modultated signal and 
the magnitude of the discrete frequency components. In this chapter, we cOB^are 
the expressions from Refs . I48jand[49j foi the spectrum of split-phase FSIC signals 
and determine the conditions under which the expressions coincide. One unique 
conclusion that can be made from all expressions is that the split-phase FSK has 
very significant discrete frequency components at the signalling frequencies, which 
are usef\a for bit synchronization. 


II . The Spectrum of Split-Phase FSK 

The paper by Hartmann [4s] has presented a general and very useful expression 
for the spectrum of split-phase CManchester coded) FSK with continuous phase. The 


main result of the paper is given by Eq. (l^) which can be siirplifted to 


^G(f) = ^ (cosirD - cos-rrF)' 


Cr^ - D^)^ 


2 CFsinTTF - DsimrD)_ T ^ n) (6.l) 




where we have followed the notations of Ref. 48 and let T be the bit duration, 

0 ) = 2lTf . , i = 1,2, be the two signalling frequencies, 2Af = ^2 “ ^1’ ^2 —^1* 

i 1 

F = (f - f )T, f =1* (f^ + ^2^’ ^ " (Af)T be the deviation ratio. The paper 

mentions the result by Shehadeh and Chiu [49] who assume that f^ and f^ are even 
integers multiplied by l/T. Under this assun?>tion, D is an integer (D ^ 0) and 
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Eq. (6.1) becomes ^ 

L„,„ . . 0.5, .(» - ») 

where p is a positive integer. 

The result in Ref. 1*9 mentioned above is the E<i. (66) which should he corrected 


to read, 2 2 ^ 

= 4 - “i’ ^ 15- “i’ I5- - “2> fs- "2> 

2 , bi \2 2 

■ (J) + ( — ) cos 

1 O 


A^t /sin^Cw^ - io)t /U'\2 
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— COS <{> 2 COs( 0 ) 2 - 0 ) )t^/2] 
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(1 + - ) (1 + r ) 


“1 “2 


( 6 . 3 ) 


where we have followed the notations of Ref. U9 with t^ T. Since f^^ and f^ are even 
integers multiplied hy 1/T, the condition is sufficient for the phase of the 
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'■ilgnal to \>e continuous. Letting ^ averaging Eq. (6.3). over <f» one 

jets after some simplifi cation , 


2 2 2 

®FSK-C^“^ = fg- «(» - “ J T?r '5(“ - + T?r " 2 ^ 


.in-(c0j_-<o)T/l»^2 ^0.^' 


-t-r (a,;-<.)T/TT j 


(1 + - f 

“l 
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. 2, 


1 + — — 

(ii-, Oi^ 
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4 (w., (oj^-a))T/4 1 2 /,.Wv/..(u 


(1 + ^){\ + - ) 


0) 


U). 


(6.M 


In the last terra, cos(wj^ - is ' always equal to one since f^ and fp are even 

integers multiplied hy l/T. Letting A = 2, folding the negative portion of the 
spectrum over to positive frequencies ^ and letting 

D = - f,)T, F = (f - f„)T, + f,) 


2^2 

one finds the nomalized spectrum, 

f°FSK-C^"> = °-5 

/■ p \ 2 

I sin 7t(F + D)/2 

'*■ V vCf +~W/2 


;in^ir(F - P)/^ \ 

ii ( F - d)72'^ J 


2 2 
(f T - D) + (f T + F) 
c c 

(2f T + F - D)^ 

0 

(f T + D)^ + (f T + F)^ 
c c 

(2f T + P + D) 

c 


2sia2.(F . D)/2 .in^w(F - D)/2 . ^ ^ 

tt(F - D ) / 2~ 


(F + D)/2 


(2f T + F)^ - 
c 


(6.5) 


The discrete portions of the spectrum are identical in Eqs . (6.2) and (6.5). Tlie 
continuous portion of the spectrum in Eq. (6.5) has heen computed hy H.P. Hartraatin 
and is shown in Fig. 19. As we can see from the spectral plots,- the two expressions 
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(Eqs. 6.2 and 6.5) agree for the case f^T » B. If this condition is not satisfied 
(dashed curves), the expression derived from Shehadeh and Chiu [U9] takes into 
account the spectral foldover, while the expression from Hartmann [^^8] does not. 

In fact, it can be shown that the two expressions coincide for f^T » B, f^T » F, 
and D being a positive integer. If one let f^T » D and f^T » F in Eq. (6.5), one 


finds, P 2 

(f T - D) + (f T + F) 
c c 

2 

(2f T + F - D) 


2(f T) 
c 

U(f T) 


= 1 
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(f T + D)^ + (f T + F)^ , (f T)^ - + (f T + F)^ 

C c _ ~ 


(2f T + F + D)‘ 


(2 f^T + F)^ - I? 
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27T 


2B + (F - B)cosit(F + B) + (F + D) costt(F - B) 


1 2 


2 2 
F - B 


( 6 , 6 ) 


in addition to 6 term. If B is an integer, Eq. (6.6) reduces to Eq. (6.2) and thus 
verifies the computer results. If B is not an integer and if f^T » B, the two 
expressions do not coincide. For example, one may consider the special case F - B. 

Eq. (6.5) becomes .22 

1 - ( sin^TTF .2 ^ 1 sin ttF x2 

“ itF (2f T)^ ^ 

c 

In Eq. (6.1), by using L’Hospital's rule, the continuous spectrum becomes 


2B^ (costtD - costtF) sinirF 


F 


= 0 




and the inconsistency is clear. 


in. Concludinp; Remarks and Acknowledgement 

We havesko^ra that the two expressions of the spectrum of the split-phase FSK 

coincide if f T »D, f^T » F and D being a positive integer. Significant discrete 
c t 

spectral components occur at; F = D, and D being a positive integer. It is assumed 

that the FSK signal has a continuous phase. 

The author would like to thank H. Peter Hartmann of AG Brown, Boveri & Cie of 
Switzerland for the fruitful exchanges on this subject which lead to Ref. 9 . 
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Chapter VII Computer Study of Phase -Locked Loop Behariors at Low Data Rate 
I . Introduction 

The performance of command and telemetry systems, useful in deep-space communi- 
cations, is frequently affected hy the radio -frequency phase error which is introduced 
at the point of reception hy means of the carrier tracking loop. In low data rate 
communications, this phase error may vary rapidly over the duration of the signaling 
interval. Causes of this type of behavior in planetary entry are turbulence, dispersion, 
attenuation and residual doppler. The phase -rariations cannot be tracked by a phase- 
locked loop of lower bandwidth, while the signal— to— noise ratio in this minimum loop 
bandwidth is too low. 

When the ratio of the system data rate to carrier tracking loop bandwidth is less 
than one, the problem of power allocation between the carrier and the data has been 
considered by Hayes and Lindsey I51j, Thomas [52]^, Sergo and Hayes [53]. For channels 
’vrith time-varying phase, Heller [5^] examined the performance' of a sequential decoding 
system. An excellent treatment of the nonlinear analysis of the phase-locked loops 
is given by Viterbi [55 J sJ^d Lindsey [56], 

In this, chapter the rapidly varying phase is characterized by a sinusoidal input 
phase, ksin (lo^t + tt/6), which models a typical phase variation in communication over 
t\irbulent media. Nonlinear analysis of the phase-locked loop behaviors in the absence 
of noise has been performed by extensive computer study of the phase-plane trajectories 
[57], [ 58 ], [59]} [60]. Readers interested in the detailed computer results should 
refer to Refs. 57-60. Only the problem formulation and the summary of the computer 
results are reported in the present chapter. Both the sinusoidal and the sawtooth 
phase detectors are considered. The sawtooth phase detector considerably disproves 
the acquisition behavior of the phase-locked loop at low data rate [ 60 ] . 
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II , The Loop Eqijation 

. Following the notations of Viterbi ( 155] , Chapter 3), ve consider first the 
differential equation of a second-order loop with perfect integrator. 






d^e. 


. AK cos<|) ^ + aAK ’sin<(> - ^ 

dt^ dt^ 

where (f»(t) is the phase error, AK is the loop gain, 6^(t) is the phase of the input 
signal, and the transfer function of the loop filter is 


(7*1) 


F(s) = 1 + - 


(7.2) 


The loop can track the frequency ramp with zero steady state error. Now we consider 
the important case that Is varying severad cycles over a bit interval of, say 

1 second which is typical in low data rate communications. The variation is normally 
caused by the time- varying channel.' Let 
e^(t) = k sin Cw^t + J) 

By normalizing the variables with 


(7.3) 


a- = ^ , X = AKt, r = 4" = y 
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m T 


Eq. (7.1) becomes 

" • ... """o • ' j. 
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in turn can be written in the state equation form as 
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(7.5) 
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where = <|>(t). It is noted from Eq. (7.5) that the larger the loop gain, or the 

loop bandwidth, the smaller the frequency of the forcing function given byEq. (7-3). 

2 

The frequency f is reduced by a factor of AK, and the amplitude reduced 

by (AK)^. In other words, the large loop gain reduces the effect of the time varying 

input phase (|)^(t). 

For a third order loop with loop-filter transfer f-uiiction 

r \ , a . b 

F(s) = 1 + - + -^ 

s 


(7.7) 


(7.8) 


the differential equation (Viterbi [55], p. 64) is 

,2 - rt 
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which jbjusing Eq, (7*4) and (7*3), can be reduced to 
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which in the state equation form becomes 
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T* 

b« = — and x_ = (|)(t), 

(AK)^ 

Phase-plane analysis of Eqs . (7.6) and (7.10) is performed by using the second- 
order Runge-Kutta method (see, e.g [59]). The computer resialts are reported in 
Refs. 57-59. A. tjq)ical computer program is given *}S Appendix to this chapter. 

If the sawtooth phase detector is used to replace the sinusoidal phase detector. 


the state equation becomes 



Xg = ~ - a' saw(x^) 

where 


koj 


(AK)' 


sin ( 


0) T 

o 

AK 



(7.11) 


and 


saw(x) = X, 1 x] ^ TT 
saw(x) = X mod. 2tt, [xj > it 

Some comparison of the two phase detectors and the tanlock phase detector has 


been given by Long and Rutledge [6l] and Uhran and Lindenlaub [62] . 


III. Summary of Computer Results 

(1) Sinusoidal Phase Detector, Second-order Loop with Perfect Integrator 

Consider first k = 0.001. The loop behavior depends on the ratio f^/AK. It is 
noted that the steady state phase error cannot be reduced to zero because of the 
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continuous input phase variation. The condition for stability is the condition for 

the loop to settle vrith a stable "limit cycle,” It is determined from the phase- 

plane study that f^/AK £ is the condition for stability. The parameter values 

for the equality to hold are called threshold values. For specified (fixed) f^, 

the threshold value of k depends on a’ = a/AK. For f^ ~ ^ ^ 

k =1 661. For fixed k, the threshold value for f is veiy sensitive to small 
th o 

parameter variation. For k = 10 , AK = 16, a = 8. The threshold value can be 

determined from Figures 20 , 21 and 22 . The trajectoiy is still converging at 

f = 0.593 . It starts to diverge at f = 0 . 59^3 and the loop is clearly unstable 
oth ’ o 

at f = 0.59^8. 
o 

(2) Sinusoidal Phase Detector, Third- 027 der Loop with Perfect Integrator 

For k = 0 . 001 , f^/AK £ 1/4 also appears to be the condition for stability. The 
increase of the parameter b' only causes the trajectory to drift more with larger 
steady state phase error [59J . A careful comparison [59] between the second-order 
and the third-order loops indicates that the third-order loop has no real advantage 
over the second— order loop at low data rate, 

(3) Sawtooth Phase Detector, Second-order Loop with Perfect IntegratOl* 

For a = 5 , AK = 16 and the specified values of f^, the threshold values of k 
can be tabulated as follows : 



k , , ( s awtooth ) 
, th 

k, , (sinewave, 
th 

1.1 

8.98 

1.661 

1.0 

10.60 

3 

0.8 

16.215 

5 

0.5 

40.5 

10 

0.25 

163.0 

50 


where the estimated threshold values of k for the sinusoidal phase detector are also 


included. 
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For input signal phase with frequency ramp of slope R rad.^ec^, Viterbi ( [55J j 
Eq. 3*27) has shown that the condition for stability for the sinusoidal phase 
detector is R ' < a* where R' = R/(AK) . By using the sawtooth phase detector, the 
condition for stability is determined from the con^uter study as R' ^ 3.017 a’ 
which is a considerable improvement over the sinewave phase detector. 
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Chapter VIII Coaing/Modiaation Selection and the Interplex System 

I . Introduction 

In this chapter vre shall be concerned vith the selection of coding/modulation 
techniques for the Pioneer-Venus main probe and small probe communication links { 63 ] . 

We then discuss a very promising modulation system, namely intei^lex, a strai^t 
PCM/PSK approach which obtains a coherent demodulator reference for binary PSK when 
the carrier is completely suppressed. 

II. Coding/Modulation Selection 

The main probe communications link, like most space applications, is constrained 
by effective radiated power but not band\fidth. Hence coding to inci*ease channel 
efficiency is desirable so long as the resulting complexity is manageable. Most 
coding techniques commonly used, however, are designed to improve efficiency of 
white Gaussian noiee channels and are extremely sensitive to burst or fading effects 
which may be anticipated in the Pioneer- Venus missions. The heavily interleaved 
block-type codes designed for bursty channels are generally not very efficient. A 
solution is to concatenate two codes; an inner convolutional code and an outer Reed- 
Solomon (RS) with an interleaving buffering scheme as discussed in Chapter V. It is 
anticipated that when the channel is well behaved (i.e. not bursty), the required 

will be 3 dB + 1 dB with these codes . During periods of deep fading most errors will 
be corrected provided that the fading characteristics have been reasonably well 
predicted and the channel is below threshold no more than 30 % of the time [6U]. 

Qgjtain convolutional codes , namely the diffuse threshold— decoded convolutional and 
Gallager's adaptive error correcting scheme are capable of correcting random errors 
and extended error bursts and are thus suited to channels with memory. To design 
such codes effectively, it is necessary to understand the fading spectrum rather 
well. Presently the precise fading characteristics are not available (see Chapter IV). 
We have to make use of the best available information in code design. 
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From the nature of the convolutional coding systems and from test results, it 
is found [65] that the Gallager code is more suitable to channels vhere bursts are 
dense and intervals between bursts are clean, and the diffuse code is optimum for ^ 
channels where bursts are more gradual and background error rate is not negligible. 

The adaptive Gallager code ’ [66 Jwhich employs convolutional coding to combat random 
errors and a form of time .diversity to combat burst errors appears to meet the needs 
for the communication situation under consideration. The reasons are: 

(1) It offers the simplicity attendant to using Viterbi decoding algorithms when 
the channel is not experiencing severe fading. 

(2) The bursty channel actually exists for a small percentage of the overall mission 
duration. The code is suitable for both bursty and additive Gaussian noise 
channels . 

(3) Test results indicate that at a burst-correction capability of about 6 seconds,, 
■thc'ccde is error-free -90- ^percent of the time which is mT.ich better than the use 

of simple diversity system. This capability is adequate for the Venus atmospheric 
channel as the average duration for each deep fade is less than 6 seconds. 

Errors in the guard space, however, must be corrected as much as possible. 
Modulation choices are not independent of the coding choice. Two types of 
modulation have been considered for use with heavily coded systems. These are 
coherent PCM/PSK/PM and noncoherent MFSK. Generally the coherent system is more 
efficient except at veiy low data rates in that, case the noncoherent system is better. 
It appears that the baseline modulation /coding scheme for the large probe should be - 
based on coherent PCM/PSK/PM modulation with convolutional inner and RS outer codes . 

Small probe communication geometry is similar to that of the large probe, 
particularly during the terminal descent phases; that the communications environment 
will be essentially the same. In particular, the possibility of deep signal fades 
due to atmospheric turbulence will be present. Consequently, the coding approach 
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suggested for the large probes, a convolutional inner code concatenated with a RS 
interleaved outer code, is also appropriate for the small probe. With regard to 
modulation, however, the low data rate for the small probe results in a poorer rat* 


of the required carrier power to data channel power, in a PSK/PM system, than was 
the case in the large probe. In this case the noncoherent MFSK system would be a 

b ett e r choi ce . 

III. The Interplex System 

It is well known that maximum efficiency of a coherent single-chennel PCM/PSK/PM 
system is achieved by completely suppressing the RT carrier, and that a coherent 
local carrier reference required to demodulate can be established by means of a 
Costas loop [67], a squaring loop [551 or another type of decision -directed tracking 
loop. However, in the existing two-channel system, it is theoretically not possible 
to completely supyress the RF carrier without simultaneously eliminating one of the 


channels. RecenLly a new two-channel modulation scheme called Interplex [68] has 
been suggested where the interraodulation loss and the RF carrier power can be 
eliminated without compromising any of the advantages of the existing system. Not 
only is the Interplex system more efficient than the existing system, when some 
unsupppressed RF power is transmitted, but it permits 100^ of the transmitted power 
to be allocated to data-bearing sidebands while preserving two-channel PCM/PSK/PM 
operation at all ratios of channel powers. In this mode, the RF carrier is completely 
suppressed; therefore, it is necessaiy to develop a method for maintaining fi*equency 


and phase sync at the receiver by methods other than the standard phase-locked loop.- 
It can be shown that although Costas or squaring loops can be used, their 
performance deteriorates rapidly as the ratio of the channel powers (or data rates) 
a = P./P. increases until no tracking is feasible when a = 1. Timor and Burman [69] 

C. J. 
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developed methods for tracking the phase of a suppressed carrier for the two-channel 
Interplex system. The resulting RF power will he independent of a. They also 
considered the performance of the system in the presence of noise. Further comparison 
between the conventional and the Interplex system is available in Ref. 70- 

The Interplex system offers a good possibility for maintaining coherent operations 
in both large probe and small probe communication links. Performance of the system 
in the presence of deep fades remains to he studied. 
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Chapt.er IX Conclijsion and Further Work 

I . Conclusion 

In his efforts to understand more about distant planets via unmanned space 
missions, man has benefited far more than the knowledge he gained about such planets. 
In the communication system design, for instance, thei'e have been numerovis problems ; 
which were never experienced before in commercial and military communication 
applications. The research and development efforts in deep space communications 
have added an important dimension to the field of communications. For the first time 
there is considerable study in designing communication receivers operating in the 
presence of large frequency uncertainty [Tl] • possibility of deep fades due to 

the Venus atmosphere has prompted careful study of the Venus channel, A need for 
high communication efficiency leads to intensive efforts of designing better codes 
such as the concatenated code, and decoding algorithms. The availability of new 
digital processing techniques has improved the software receiver design. These 
have been the topics in the present report. 

It is reasonable to conclude that, although there are still many problems to 
be studied, we now have a much better understanding of the. low data rate digital 
space communications. Continued efforts in these studies will assure the success 
of the forthcoming tintried deep space missions. And many of the results we obtain : 
from these studies will, also be very useful for solving the communication problems 
on earth . 

II . Recommendations for Further Work ^ 

In Chapter II, the software or hardware inplementation of the proposed all- 

digital system is recommended. 

In Chapter III, re-evaluation of the system performance is recommended when 
more precise fading characteristics become available. 





In Chapter IV, precise fading characteristics remain to be examined. 

In Chapter V, the best available performance from using the concatenated 
code (Code III) must be eveiluated theoretically. 

In Chapter VI, further study is needed on the power spectral density of the 
angle modulated signal with split-phase baseband format. 

In Chapter VII, theoretical analysis of the phase-locked loop behaviors with 
rapidly varying phase error is much needed. Further computer study on the third- 
order loop is also recommended. 

In Chapter VII, the modulation selection for small probe communication link 
needs further study in view of the available Interplex system. 
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0. 29329 
0.20870 
0. 17916 
0. 1 3962 

■ 0. 10508 
0.07059 
0.03600 
0.00196 

-0.03308 
-0.06762 
-0. 10216 
-0. 13670 
-0. 1/129 
-0.20578 
-0. 29032 
-0.279B6 
-0.30990 
-0. 39399 
-0. 37898 
-0,91302 
— 0. 99756 
-0.96210 
■-0. 51669 
-0.55116 
-0.58572 
- 0.62026 
-0.65980 
-0. 68939 
-0.72388 
-0.7589? 
-0.79296 
-0.82750 

■ ■ - 0.2 


IPg.il 


piiAs. " P/tfifie plot i 

£jrC7.6) 1- 

~~ i ^~ rT 6 [' 

f, = o.5-?43 1 





- 0 . 2 ' 
2.9613^ 
2 . y 8 '+ 3 5 
2.80736 

■ 2.73037 
2.65338 
2. 57639 
2.A99A0 
2 . 422^1 
2 . 3 ^ 54 ? 
2.26343 
2. 19L44 

■ 2. 1 1445 
2.03746 
1.96047 
1.88348 
1. 8'0649 
1.72950 

■ 1.65251 
1. 57 5 52 
1 . 49853 
1.42154 
1. 34455 
1.26756 

■~1. 19057 
1. 11358 
1. 03659 
0.95960 
0,88261 
0,80562 
‘ 0.72863 
0.65 164 
0.57465 
0.49766 
0,42067 
0. 34368 
0. 26669 
0. 18970 
0. 1 1271 
0.0357? 
-0.04127 
-0. 11826 
-0. 19525 
-0.27224 
-0.34923 
-0.42622 
-0.50321 
-0.58020 
-0.65719 
-0.73418 
-0.81117 
-0.888 16 
- 0.2 


t 
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MA IN PGM 


IMPL rClf REALMS 
DIMENSION A ( 1.00 
READ( 1 , 10 
FORMA n’ 2 \ '5 f F5.0' 
MB-F-I 

2 ❖^MBITS'"' 

NB~N-1 

W RITE ( :5 t 2 it ) M » N « i 
FORMAT ( ® SIT 
resvs) psisx' 
IM 3 ) = I «-P { l_) 

AL9 = 0. 

IF (M-i) WCS 1 70, ' 
DO SO L=^2,M8“ 
XL==L 

AL9-^.^L9->''ULaG(XL 
DO 200 I==Lr3 

po-p { I J 

C2 = P9^'=^ { I . /NB ) 

2 1-0 ■ 

CALL GAHlNlMtMB 
Fl='( H-^C2 ) /DEXP ( i 
Z l-Zl-El 
I FI DABS { ET ) .GE ; 
XC I ) -Zl 

WRITE(3,60U>m' 
FORMA r I IXf IP2DI 
HR 1 1 E ( .3 , 70 ) 
FORMAT ( * OF /NO * , 
K6-NBI I L + l 

09- 1 . +GMSQ+R6 
0 9- ( L T+GMSO J /09 
A2-M-xR6^'^GhSQ/D9 
HI- ( XI 3)-X{ 1 ) r/l 
N 6 L - N 7> +• 1 
DO CO-H) ’ I-l»06 r“ 
X9-XI i )+HI^'{ [-1 
C.ALL OSUBMl Mr A.2~ 
CALL GAMI.NIM,MB 

H-H 8 - { ( i\io - L V 

n - n 1 / 1 .: a r v rl o u l w 

A(I)-Q^<H .' 

AIl-0. 

UG S /O' I-2VN6V2 
AU-A I S A ( I ) . 

A I 1-2 I 1 + . ( 

N6M--06-I 
UO 710 r-3',N6'MX* 
A 1 l-AI H-A ( [ ) 


-HrO-Z) 

3 1 1001 ) /P I S T XI 3) 
»)M,NBI rS_rGMSGrPt Dp 
1 1 0 • 0 t I 3 .3 F 1 0 » 0 ? i j 1 0 • 


0 r N 6 

r * N- S i 7 , SX r * GM 
,{/)■ ‘INITI AL*7)'" 


P ( 2 } r 
0) 


E6, E7, E8 , ERROR 


, F 1 0 . 


oxD wr 


IS,» QRDINA 


. t H I 

M3Li)GTZD~Z1-AL9' 
I-63GQ ■■T0‘L8C3V'"'"‘" 


D PE ( N6/2 ) S lOX 
,aT7('M>i^GMSQ) 


l.Xt *PES lAX, 'ERROR 


’=iPB'9'rERR0R',‘in' 



0 4 V r ~4 A 1 o ^ 

, *7 ^ A -7 — U 7 / 


3‘+'A ( N6 D 




30 7 


SUBRUUT ItviC- QSU3M ( M.9» A2 t B? r 
COMf^’UTES Tht: Gh^jeRALl^C-0 0 FUNCTIDN 
M9- M t A2" A-:";=?./2 * B.2= B;^*2/2 p 

I Mf>L 1 C I T Kt AL-:' 3 { A-ii fV-L) 

XKO-o\iXP { -\2) 
l)O = 0l:XP (-B2 ) 

>0=L,*-QO 

I F ( r 7 »GT • I ) GU TO 316 
^1^.5^'^SQRT { 2,-A2^'!3_2 ) 

■ Q=U-GO^XKQ 
K — I ♦ 

■ XK0-A2^:^XK0/R “ 

nO = H2"^DO/i< . __ 

CO^GO-DO 

6A = G0-XK0 ■ __ _ . 

C-Q-F^t 

I r- ( R . L T . R I ) GO T O 314 

I F ( 04 • L T . E 9 > GO 'T O' 3.3 '{ 


•^HERF 

E9 = ERROR 



3 1 4 


316 


10 


11 


R — R + 1 — 

GO ]0 30 7 

0B = 09-1 

XOO-G3 ■ 

{ I . +XMH=^^( DSQRT ( 1 o A2 -3 2/ l_XMc^^2 

0 3-00 

no 10 M7=LfH8 , „ 

0 3==b2-D.3/M7 

GO^GO-D3 
iJlJ M H=^2. tM9 

00 = 00^:^32/ { M~ U _ 

(j = U-GO^:'XKO 


))-!.)) 


R=^l * . . 

320 XK0 = A2-XK0/R‘ 

D0=B2-D0/ i R+XMR ) 

G0=G0-00 

E^xGO-XKO 

fV' ( R. LT • R I ) GO rO 3.35 
1 F( E4.LT'. E9VC.0 rO"'33 /' 
335 R^R+1. 

GO ro 32 B" 

337 KETURM _ __ 

■ ' END ■ ' ■ 


GAMIO_^ 

.subroutine GAFliiM(M9tNB,X9,H5) _ 

B”' CORPUTbS NURRiALI /ED INCOMPLETE GAMMA FUNCTlDH WHERE 

^ . M9 r: HB=H~1 D X9 =X « H5 = THE ANSWER 

[MPlICIT REAL=:'a’( A-HtU-Z L. 

n3 = DEXP ( *- X9 ) 

H5=l.-05 

{ F ( M 9 4 E Q o 1 ) GO TO ’28 7 • . 

00 286 

05 - ( X9/M ) ^!>5 
2B6 H5-H6-D5 

287 return 
ENI.> 


A2 

= THl; ANSWER 


/13 


..V I l=^z, . 

UO 7 70 I'-* I, NO 1,2 
770 n( [/2 + U^A( 1 ) 

A I 1^0. . ^ 

MO- No/ 2 

DO B20 1-2,N5t2 

8 20 Ari = AU^0ir) 

A f I = 2 . A 1 L +•. SV( B f 1 ■■) ^0 ( N 5 + i ) )' 

N5M-M3-1 , _ _ 

DO 8 60 I-J,N0m,2 

860 AII=A11+D(I) 

A t i-A. ): 1/3. 

A [8- ( A 12“ A I 1 ) /I 3. 

A1 3-A 12-t'A [ 8 ‘ ” 

A I I = 1 . - A 1 1 

A 12-1 . -A I 2 

AI3=i.-AI3 
A!'B = -AI8 

WK I T E ( 3 , 9 0 ) F6 , A I I , A 1 2 , A I 3 , A I B 
90 FODMAT ( IX , 1P0018. 6 ) 

.E iJ = E b + E B 

TF(E6“E7)300,'300, 1 

3000 CALL EXIT 

END 


APPENDIX - CHAPTER 


MAINPGM 

llTLi: : ANPLITUCE FADING STATrSTICS AFS-l. 


C 

%V a, ^ vTf N»^ nV ->f «^V 

'V^ ■'«'• T- •%- V 't" •'}' '•*' N' 

: c 

IHIS FPGGRAN 

uC 

PkCBAO IL I TIES 

c 

FCP EACH FHCB 

c 

LEVELS IN OB 

c 

«R2’ ARE THE 

c 

FERENT SC ALIN 

c 


c 

IN THIS IM^CGR 

c 

FORMULARS BY 

c 


c 

FREQUENCY - 2 

c 



ijc V V **v A* ^ V , V 


R F A ! K 1 - K ? - K 


FGRNAT ( F20 , 3 ) 

2 

FORMAT ( GX » F12 

3 

FORMATI’ P-', 

4 

FORMAT ( lOX, « Z 

5 

FORMAT I ’ 0 ’ ) 

10 

READ! 1 , 1 ,ENO= 


WRne( 3,5) 

PO 

torn F { 3 , 3 } P 


IvR I TE ( 3 , A ) 

30 

CALL SUBKYY, 

40 

P= I ,“P 

^0 

00 190 1=1,51 


Z=I-I, 

6 0 

H=l./*0785 

' 70 

E=EXP ( -Z/M) 

1 80 

Sl=2c 87-^E 

""90 

S2= U 95=!9' 

100 

Kl-SQRT ( EXP ( 2 

rn'o 

K2 = Sv;rT( EXP{ 2 

120 

K = Kl 

"l30 

CALL SUB2(R,K 

» 140 

RI=R 

c; n 

K = K 2 

160 

CALL SUB2(R,K 

^'170' 

R2 = R 


VO U T E { 3 , 2 ) Z , 


CON ri NUE 

200 

GO TO 10 

300 

CALL EXIT 


END 


!\A' I 

i 

\ 

I 

. \ 

I 

t 

I 

i 


i 

* I 



SUBPtCUTINE SUril{YY,P} __ 

C ij! 5.V 4: 4; ?5t li: s;< 4: «: sft jJ: >I< <? r)i ❖ s!< 'X s}f 5^ ❖ ^ -I' ❖ v >i'' ^ 5r >;« ❖ ❖ sji ^‘i i,\ j|« i> ^ sj! ^ :❖ 5j? ^ >5^ ❖ ❖ ’r 4^ >1' <' V- 

c 

C^ SUdKGuriNfc TO FHmD Y, k-^'HERE Y IS THE LEVEL SUCH TliAT FOR A GAU5- 
9 CR t\CRHAL P.ROBABlLirY D I S T R 1 5U T I CN » P=PRG0(R’<Y) 

C ’ 

C i: 4 4 4 4 4= 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 * >:? 4 * 4 4 4 4 4 4 4 4 4 4 4 4 4 -I; 4 4 4 '4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

970 C=l-“P 

9bo if{ c.u r, .:» } so rc io;>o 

990 A^l. 

UK)0 T^SQRT ( ALGG{ U/QY*2 ) ) ' . 

101.0‘ Y = 2*S!. SS 17 + T'4 ( 0^802 89 3 + 0.0 10 32 ) 

__i020 Y = Y/1 K + T4( I, ^32788 + 14 (C. 189269 + 0.C0L308^-'T1} ) . 

1030 Y-Ari T-Y) ■ 

_YY = Y _ _ 

'’l040 GO in 108 0 ' 

i050_ A^7l. . . - . 

1060 (0=^L*“-Q 

10_70 _00 TU 10 00 ^ 

■'TobO" RrfURN 

END 

SUHRUU r IJnIE _SUB2 ( R tKt YY ,P ) .. 

44T;;4 44 4^4 444444 444 44 4444 444 44 4 44 4 4 44 4 4 44 444444 4444444 44 4 444444444 4 4 444 4 444 4 4 44 

SUBRnUT I«YE' TO C'^^LCULArE T l^E 'f'a 0 1 NG LEVEL IH D6 COMPARED TO THE 

F^REG SPACE LEVEL. THE. CAL CULATJ CR^J S_8A SED.. GN TH E, MGR TCi\ ‘ S EQUA-, 

riciN'S. 


; 4 i,i 4 4 4 <x jj: 4 4 4 4 4 4 4 4 4 4 4 4 4 44 4 4 4 4 4 4 A- 4 4 4 4 4 4 4 4 4 


: 4 4 4 4 


^ IP/. ./W -j. *^n 


REAL K , K 
LOG! AA ) = 
V :^Y_ 

I r V K . G T o 
K8=K442 

'I fIk^gk 

R=K4Y/1 . 
R -8. 6858 
R=R-3.68 
"go TO 9 A 
R~ I . +K8/ 
R=R+V4K/ 
_R = R + Y442 
R=R-K448 
8-8.6858 
R5-8 .685 
R-'^R5 
GO ru 9 A 
X — LOG( P 
TF( K.GT. 
K 8-1 . / ( K 
R- 1. + K8^:' 

R - R IL!^ H 4 
■r = K'44 24X 
R=A. 3A29 
1R5-8 .685 
R=R-R5 
CrO TO 9 A 
R -8 .6853 
'R = R+Ao 3 A 
R5=n ,685 
R=R-R5 
RETURM 


8 T L 

alogVaa') 

r/) GO TO 

TlT‘"‘G(T"tQ” 

A1 A4 ( 1 , -3 
94 ( R- (K4Y 
5894LCG (S 
0 

A_._4 U . ^ K 8 
I . A I A 4 ” ( I . ■ 
4KB 442/2 A 
4 V443/ 32 . 
94L0G ( R ) 

8 94LCG (SC 


720 

. /B . 4K8 ) +K8/ A .4 ( 1 .-K8 
)4424 ( 1. /A .“K4Y/6./I . 
CRT ( l_.+K8 ) ) . _ , 

/2A.4J I. 3 54K 8 ) } 

- K S 7b 4 { ‘ L /-■ 3 . / 16 , "4 K 8 ) 
.4{1,-1Y,/A0.4K8 ) 

/ i .A LA + K'’'4 64Y44A- / 80* 

RTU.+KB) ) 


/ 6 . ) 
A 1 4 ) 


10. ) GO T 
442 ) 

[ 1.+.54K8 

4A )_/2A . 4J 

4i7' 

4{loG(R) ) 
894LCG ( SQ 


94L0G ( K ) + 
29/ ( K442 ) 
B94LCG ( SC 


0 900 

4{ 1.-X/2. )+k'B 442/6.4{" 
l.-17./2.4X+Ag,/6.4X4 


RT I i .+K442 ) 1 


A.3A294U)G(X) 

4{ I . -X/^i . / ( K4 4"2 ) ) 
RT ( J . +K4 42 ) } 


1.-2 

4 2-3 


. 54X + 2 
3 . / 2 A . 


. /3. 

4 X 4 4 


4X44; 

3 ) 
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M,\ l/jPGM 


AQ 



PL A. ‘IP 


A PAL VS IS 


v;irr( PPKPCCV 


i^C 

F I i. T 



•H-AL K 

01 Vl'-;S I0:i Y ( 1.000) , Z ( I 000 ) 
0 iP‘-"iS i.n,N FU(2) 

DATA Fn/ . ‘>04, , ‘39A^/ 


, TAU I ( 1000 ).T OTAUl ( 1000 ) 


v%‘ = 2 . ( \ . 

O^H/ ( 2 1 ■ 

P 1 = 3 , 1 . A I 092 65 

FUR(V.\r ( 4F5 . ?- ) ' ’ ■“ ■ 

FORK AI { IHO , 3X , ’ liUMBFO * » 2 1 X , * X 1 ' , 24X , * X2 ’ , 
FORMA r ( 5X , I 4 , 17 X , E 14 » 7 , 13X , E 14,7 , 12X , E 14. 
FORMAT ( IHl > 

FURP,\ I ( IHO) 

PvEAO { 1 , I , L M D= I 0 ) X I A , X 2 A , A , 0 

no 9 J = I r 2 . - - - 

K = 1 0 , 

AH=0/A ■ ■ " ' 

FRI:C = FO ( J ) 

WZ=2. I-FRtC 


2 2X , ' OTAU* , 2 0X , » TAU< ) 
7,9X,E14;7) ' 


X 1 = X I A 
X2=X2A 
TAU=0. 
DTAU=0. 
X U)=X2 


) X 2 “ A P ST ITT X-1- ) - ( K ^ ( W I 


20 / ( 2 ) ■} * S I NT ( V>) Z - T AU / A 1 ^ P I / 6 - ' 


E ) 

on H 1 = 1, 1000 ■ ‘ 

JF( I ,LE. 2 ) so TO 7 

DTAU= \13S { { Y ( I “D'-Y ( I -2 ) ) / ( (‘ZTT - I ) +Z ( 1-2 ) ) / 2 . ) ) 


TAU= rAU + DTAO 
XlH=x L4-R^.'XID 


X2H = X2 4-R- X2D 

X10H = X2H • - ■ 

X 20H=~C OS ( X IH ) X2H-A S 1 M ( X IH ) 

SI /6. ) - - 

X 1 = X 1+ H-M ( I . - 1 0 + X I OH ) 

X 2 = X 2 4- H"' ( ( 1 . - W ) - X 2 D + 04 X 2 OH ) ' " 


( K =!= ( W Z’- - 2 ) / ( A* -•!' 2 ) ) =5^ S I N { ( W Z -T AU / A ) +p 


X1D=X2 

X 20 = -C 0 S ( X 1 ) X 2- A P S I N ( X 1 ) - ( K ( W I 


2 ) / { >1= 2 ) )M< S I N P( W AU / A ) + P I / 6 . 


S) 

■Y( I )=X1 

Z ( 11-X2 

TAUl { n = TAU •“ •” 

OTAUl ( I ) =OTAU 

CONTINUE 

CALL PLOMT (Y,Z, 1000) 

vJRin:(3t‘^) " 

Hil ITE (3,2) 

CONTINUE 

CALL EXIT 

lMO ' ‘ ' 


